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ABSTRACT 

In magnetohydrodynamics (MHD), the magnetic field is evolved by the induction equation 
and coupled to the gas dynamics by the Lorentz force. We perform numerical smoothed par- 
ticle magnetohydrodynamics (Spmhd) simulations and study the influence of a numerical 
magnetic divergence. For instabilities arising from V • B related e rrors, we find the hyper- 
bolic/parabolic cleaning scheme suggested by iDedner et al.l (l2002h to give good results and 
prevent numerical artifacts from growing. Additionally, we demonstrate that certain current 
Spmhd implementations of magnetic field regularizations give rise to unphysical instabilities 
in long-time simulations. We also find this effect when employing Euler potentials (diver- 
genceless by definition), which are not able to follow the winding-up process of magnetic field 
lines properly. Furthermore, we present cosmological simulations of galaxy cluster formation 
at extremely high resolution including the evolution of magnetic fields. We show synthetic 
Faraday rotation maps and derive structure functions to compare them with observations. 
Comparing all the simulations with and without divergence cleaning, we are able to confirm 
the results of previous simulations performed with the standard implementation of MHD in 
Spmhd at normal resolution. However, at extremely high resolution, a cleaning scheme is 
needed to prevent the growth of numerical V • B errors at small scales. 
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1 INTRODUCTION 

On large scales magnetic fi elds are obse rved within galaxies, along 
galactic outflows (see e.g. lBeckll2009l. for a review), within the 
intra-cluster medium (ICM) (iGovoni l2006h and ev en in the fila- 
mentary structure of the Universe ( Dolag et aDl201lh . Faraday Ro- 
tation Measures (RM) towards extended radio sources in and be- 
hind galaxy clusters (used to probe the magnetic field within the 
ICM) can demons trate the turbu lent nature of the magnetic field 
in galaxy clusters (lTribblelll99l1) . Such observations, yields typi- 
cally /iG values for the central magnetic field within galaxy clus- 
ters, whereas the reconstruction of the magnetic field mainly leads 
to a radia l decline, consistent with a turbulent magnetic field dis- 
tribution teuchar & E nBlinll201ll: iGovoni et aLll201Ql : IVacca et all 
[2012; Guidetti et al. 201Ql: lBonafede et aLkOlof) . 

The presence of turbulence within the ICM is expected 
from hydrodynamical cluster simulations in a cosmological con- 
text, predi cting a significant amo u nt of turbulence within the 
ICM (e. g. 'Brvan & Norman "l998'; 'inogamov & Sunvaev '2003|; 
Rasia et a l. 2004; Dolag et al. 2005; Vazza et al. 2006, 2009, 201 1]; 
lapichino & Niemevejl2008l : iPaul et al.ll201ll) . Therefore, the tur- 



bulence should amplify the magneti c field within the cluster a tmo- 
sphere towards equipartition values (ISubramanian et al 1l200g) . 

Non-radiative cosmological simulations of galaxy clusters 
following the evolution of a primordial magnetic seed field 
have been pe rformed using smoothed particle hydrodynamics 
Spmhd codes JPolag et al.lll999l . l2002l : IPolag & Stasvszvnll2009l : 
iBonafe de et al."2011) as well as adaptive mesh refineme nt (AMR) 
codes (Bruggen et al. 2005; D ubois & Teyssier 2008; Co llins et"an 
|20i3; IXu eTall |2010|; lMiniati & Martini l201lh. Such cosmolog- 
ical ideal MHD simulations also predict a n iagnetic spectrum , 
efficiently redistributing the intern al energy (iDolag et al.l l2002l : 
iBruggen et al.ll2005l : IXu et"aDl2009h into the ICM. 

However, to capture the features characterizing the magnetic 
field within galaxy clusters, their build-up process within the large- 
scale structure of the Universe has to be followed correctly. There- 
fore, numerical codes have to resolve many orders of magnitude in 
dynamical range. In this process, a variety of dynamo theories will 
be involved and tested, whereas it can point out towards a turbulent 
dynamo or shock enhancement of magnetic fields. 

Here we present a n extension of the standar d Spmhd imple- 
mentation in Gadget (iDolag & Stasvszvnl 
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is on simulations without radiative processes, to gain deep insight 
in the performance of the code itself. Additional physical processes. 
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would change the behaviour of the system and make it difficult to 
disentangle the possible small-scale effects resulting from a non- 
vanishing V • B the case of our ideal MHD simulations. With this 
cleaning scheme we are able to perform extremely high resolution 
simulations of galaxy clusters in a cosmological context. By reach- 
ing a spatial resolution of kiloparsec in the center of the cluster we 
are able to directly compare the structure of the magnetic field ob- 
tained in the simulations with the observed structure on the smallest 
scales observed. 

The paper is organized as follows. In section |2] we present the 
further development of the MHD Gadget code. Results from typ- 
ical standard tests are shown in section |3] In section]?] we present 
simulations of galaxy clusters and compare the predicted magnetic 
field structure to observations, before we summarize and conclude 
in section [5] Additionally, in the appendixes, we show more de- 
tailed the effect of dissipation and resolution on the results. 



When reading output data from the simulation the magnetic field is 
converted back into physical units by multiplying it with l/o? . 

To capture in more detail situations where the magnetic field 
structure is folded on the resolution scale, an additional time step 
criterion for every particle i can be constructed: 



hi 



typical 



hi 



47r(V • B) 



(2) 



where h is the Spmhd smoothing length. This criterion allows to 
capture regions of high numerical divergence, where the magnetic 
field structures reach the resolution limit and put the particles on 
lower time- step ensuring a more detailed evolution of the magnetic 
field. Although this additional criterion rarely overcomes the stan- 
dard time- step criterion, there are situations where it seems to be 
quite helpful to follow the local dynamics with more details. 



2 SPMHD IMPLEMENTATION 

We are starting from the Spmhd implementation of 

jPolag & Stasvszvnl l2009h within t he cosmologi cal N- 

Body T r eePM/ Spmhd code Gadget ISpringel et aD (l200lh : 
Springell (2005h. For the details of the implementation we to 
Dol ag & Stas vszv or to a more general review on Spmhd 

by [price ( 20 ij^ 

In short, the evolution of the magnetic field is directly fol- 
lowed with the induction equation. The magnetic field acts on the 
gas via the Lorentz force, written in a symmetric conservative form 
using the magnetic stress tensor. Also, as originally suggested by 
IPrice & Monaghan ( 2004) the fast magnetosonic wave velocity re- 
places the sound velocity within the computation of the signal ve- 
locity controlling the artificial viscosity and the time-step. In the 
calculation of the gradients and divergence estimators, we follow 
the standard Spmhd implementation. 



2.1 Instability correction 

To take into account the tensile instability in Spmhd, which oc- 
curs when the magnetic pressure exceeds the gas pressure and the 
force between particles is becoming attractive, a correc tion term in 
the fo rce equation is used. This t erm - introduced by lB0rve et aP 
(l200lh and further developed in lB0rve et al.l (l2006l) - subtracts 
from the equation of motion any unphysical force resulting from a 
non- vanishing numerical V • B. Contrary to the original implemen- 
tation, we restrict the correction to not exceed the Lorentz force, 
which is necessary at strong shocks. Therefore, we evaluate the 
correction contribution and if necessary renormalise to be only as 
much as the Lorentz force. 



2.2 Time integration 

In lDolag & Stasvszvnl ( [200^ the evolution of the magnetic field B 
is done in physical units, so that in cosmological simulations the 
induction equation contained a term of —2B to capture the cosmo- 
logical dilution due to the expansion of space. However, defining 
the magnetic field in comoving units {Be = B/a^), with a the 
cosmological scale-factor allows to drop the —2B term and the in- 
duction equation becomes 



dBc. 
dt 



2.3 Divergence Cleaning 

It is of fundamental interest in Spmhd simulations to keep the 
magnetic divergence V • B arising from the numerical integra- 
tion schemes to a minimum. We implement in to the GADGET code , 
the divergence cleaning scheme introduced bv lPedner et al.l(l2002h . 
which evolves an additional scalar potential representing non- 
vanishing V • B introduced artifacts. By construction, this potential 
^ propagates the numerical errors outwards the simulation, while 
damping them, by subtracting the gradient of in the induction 
equation. This method is also widely used in Eulerian codes (i.e. 
I Cecere et al.||2008|; I Anderson etaPbOOd: iKeppens et aLll2012h and 
recen t ly introduced in th e moving mesh code Arepo ( Springell 
l20ld : IPakmor et al.ll201lh . First attempts t o use t his technique in 
Spmhd were made bv IPrice & MonaghanI (l2005h , being not sat- 
isfactory. They found only a mild improvement of the numerical 
V • B errors and in some test cases the cleaning scheme was even 
causing instabilities particularly in 3D. The implementation shown 
here, lowers the V • B error in general, does not show unwanted 
effects and addi tionally leads to a very small numerical diffusion. 

Following IPedner et al.l (|^02), therefore assuming a non- 
vanishing V • B, an additional term entering the induction equation 
can be derived: 



dB 

~dt 



(3) 



To be energy conserving, the removed of magnetic energy is trans- 
ferred into internal energy or entropy A at a rate of 



dA 



7- 



B^ • (VV^)i 



(4) 



where (7 — l)/p^~^ is the conversion factor from internal energy 
to entropy with the adiabatic index 7. 

However, the scalar potential iJj has to be chosen in a way to 
actually removes numerical errors, bedner et al.l (I2OO2I) found that 
the most effective solution is to construct and evolve i/j propagat- 
ing the errors away from the source (i.e hyperbolic cleaning) and 
damping them (i.e parabolic cleaning). This results in the following 
evolution equation for i/j: 



dipi 



(ch)?(V-B).-^ 

Ti 



(5) 



(Be- V)i;-Be(V' 



(1) 



which shows that 1^ now satisfies a wave equation propagating the 
errors outwards the source with a speed of Ch (first term of the 
equation) and decaying them on a timescale of r (second term in 
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the equation). It is again natural in Spmhd simulations to relate 
the propagation speed to the fast magnetosonic wave, hence using 
Ch = crv™^^. Note, that this velocity does not have to be related 
with any special quantity per se, even can be a constant value. Also 
the timescale r can be related to a typical length scale (smoothing 
length h) and velocity resulting in h/Xv^^^, only leaving dimen- 
sionless numerical constants A (parabolic), a (hyperbolic) of order 
unity. We choose values of A = 4 a nd a = 1 to recover the best 
solution in lPrice & MonaghanI (l2005h . 

Similar to the tensile instability correction of the Lorentz force 
(lB0rve et al ]|200lh in the equation of motion, this method can lead 
to instabilities. In particular, situations where small scale structures 
in the magnetic field lead to an imprecise calculation of the V • B 
source term. This situation will manifest as an over-correction of 
the induction equation. Hence, we use a limiter for the cleaning 
contribution in a similar form, not allowing the correction in the 
induction equation to be larger than a given value Q weighted by 
the local induction value. When the correction exceeds the original 
term, we renormalise it. To ensure stability, this ratio Q has to be 
less or at most equal to 1. Testing different parameters we found 
a value of Q = 0.5 to be sufficient in ensuring a proper evolution 
of the magnetic field while avoiding over-corrections due to the 
cleaning scheme^ 

Note, that lPakmor et al.l (l201lh use a global maximum for the 
signal velocity in the evolution of the scalar field V^. However, our 
Spmhd formulation of the cleaning equations takes all quantities 
per particle, in order to avoid any dissipation or overcorrection in 
already stable regions. 



3 APPLICATION TO TEST PROBLEMS 

The problems in structure formation are very complex and astro- 
physical objects of interest evolve in a strongly, non-linear way 
from the initial conditions to the final stages during the different 
cosmological epochs. To be confident about the numerical results, 
the hydrodynamical solving scheme has to be tested properly and 
compared with known analytical solutions. Therefore, we tested the 
new cleaning scheme in an extensive series of shock t ubes and pla- 
nar tests, similar as done in iDolag & Stasvszvnl(l2009h . In the same 
way as done previously, we performed all the tests by setting up a 
fully three-dimensional glass like particle distribution, to obtain re- 
sults under most realistic possible circum stances and c ompare then 
with the solution obtained with ATHENA dStone et al.lE008.) (in ID 
or 2D, respectively). Additionally, the V • B errors are defined by 
the dimensionless quantity 

ErrvB, = |V-B.|-^ (6) 
\tii\ 

which can be calculated for each particle, and determines the relia- 
bility of the results. 



3.1 Shock-tube tests 

The most commo n MHD shock tube test is the lBrio & Wul (Il988h 
shock (Test 5A in IRvu& Joneslll995h . In this shock-tube, a shock 
and a related rarefaction are moving togethe r. For the standard 
scheme, we used the method described in iDolag & Stasvszvnl 
(l2009h . Additionally, the instability correction as described in sec- 
tion |2] was used. The new Spmhd results generally agree with the 
solution obtained with ATHENA, although there is some residual 
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Figure 1. Shocktube 5A at time t = 7 with the Dedner cleaning scheme. 
In the first row is the density (left panel), total energy and pressure (middle 
panel) and the x velocity component (right panel). The second row shows 
the y velocity component (left panel), the y-component of the magnetic 
field (middle panel) and the measure of the V • B error (right panel). There 
we also show the V • B error of the Standard MHD implementation (in 
blue). The colored lines with error bars (blue for the pressure in the upper 
middle panel) show the Spmhd results, the black lines are reference results 
from Athena in a ID setup. 



scatter in the individual particle values within the 3D volume el- 
ements, as well as some small scale noise, especially in the low- 
density part. Note, that the mean values for the internal energy, the 
velocity or the magnetic field, can locally show some minor but sys- 
tematic deviations from the ideal solution. However, the total en- 
ergy shows much better, nearly unbiased, behaviour. This demon- 
strates the conservative nature of the symmetric formulations in 
Spmhd. 

As can be seen in Fig. ([T]), the cleaning scheme achieves a vis- 
ible reduction of the V • B error and lowers the numerical noise. 
To find good numerical parameters for the cleaning scheme, we 
performed a parameter study test series. This study takes into ac- 
count the estimated noise and accuracy, as previously explained in 
iDolag & Stasvszvnl (l2009l) . The best results are obtained for pa- 
rameters ranging from A = 2^5tocr^l. We therefore 
have chosen the va l ues of A = 4 and a = 1, similar as used by 
IPrice & MonaghanI (l2005h . These numbers imply (see Section |23] ) 
that the V • B errors propagate at the signal velocity and will be 
damped within approximately 2^5 smoothing lengths distances. 

Fig. ^ shows a case, where the original cleaning scheme 
without limiter fails. However, with the limiter Q ^ 0.5, no insta- 
bility can be seen and we notice the V • B error again to be smaller 
than with the Standard implementation. Comp aring results of dif- 
ferent regularization schemes as presented in iDolag & Stasvszvnl 
(|2009f) and the cleaning scheme presented here, we note that the 
cleaning scheme does not smooth sharp structures within the dif- 
fe rent shock-tube test s , as o ther regularization methods described 
in lDolag & Stasvszvnl (l2009l) are doing. 

To check the performance of the new implementation in de- 
tail, we ran the full set of different shock-tube tests as presented in 
IRvu& Jonesl ( ll995h . The results are summarized in Fig. (|9), where 
the black diamonds represent the results of the 1 1 shock-tube test. 
The cleaning scheme is stable for all 1 1 shock-tube tests and the 
V • B values decrease on average by a factor 2 3. 
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Figure 2. Shock-tube IB from lRvu & JonesI ( 1 19951) . In black are the result 
from Athena, in blue our cleaning scheme with Q delimiter and in red 
without. From upper left to bottom right, density, x velocity, y magnetic 
field and V • B errors are plotted. Again, in the last panel, are V • B errors 
for the Standard case. With the limiter Q ^ 0.5 the cleaning scheme leads 
to less V • B error than the standard MHD implementation (black line in 
the lower right panel). 



V • B error by a factor of two, as can be seen from the according 
data point in Fig. 



3.2.2 Strong Blast 

The Strong Blast test is of an explosion of a circular hot gas blob 
within a static magnetized medium and i s also commonly used 
for MHD code v alidation (see for example [Londrillo & Del Zannal 
l200Ql :l Balsara & Spicer 1999). It consist of a background medium 
with a constant density p = 1, where a hot disk of radius ro = 
0.125 and pressure Pd = 100 is embedded in gas with a pressure 
of Po = 1. Additionally, there is initially an homogeneous mag- 
netic field in the x-direction, with a strength of Bx — 10. The 
system is evolved until t = 0.02 and an outgoing shock wave de- 
velops, travelling not circular, but along the magnetic field lines. 
Fig. ^ shows the density at the final time, comparing the ATHENA 
results with the Dedner scheme in GADGET. There is no visible 
difference between the Spmhd implementations and the Athena 
results. Besides some very small variations, there is no significant 
difference between the various Spmhd schemes, and all features 
are well reproduced. In this test problem, the V • B cleaning scheme 
has effectively a similar dissipative behavior than the scheme with 
artificial dissipation. As can be seen from the according data point 
in Fig. 12]), the cleaning as well as the other regularization methods 
give moderate improvements in the numerical V • B errors. 



3.2 Planar Tests 

Besides the Shock- tube tests described in the previous section, two 
dimensional (e.g. planar) test problems are an excellent test-bed 
to check the performance of MHD implementations. Such higher 
dimensional tests include additional interactions between different 
evolving components with non-trivial solutions. They can be quite 
complex (with several classes of waves propagating in several di- 
rections) such as the Orszag-Tang Vortex or simple (but with strong 
MHD discontinuities) like the Strong Blast or Fast Rotor. 



3.2.1 Fast Rotor 

The F ast Rotor test was first introduced by iBalsara & Spicerl 
(Il999h . to study star formation scenarios, in particular the 
strong torsional Alfven waves and is also commonly used 

to val idate MHD implementation s (for example see iTothI 

200QI: iLondrillo & Del Zai^ l200d: I Price & MonaghanI l2005l : 
B0rve et al.ll2006l : iDolag & Stasvszvnll2009h . The test consists of 
a fast rotating dense disk embedded in a low density, static and uni- 
form medium, with an initial constant magnetic field along the x- 
direction (e.g. Bx — 2.57r~^/^). In the initial set-up, a disk with ra- 
dius r — 0.1, density p = 10 and pressure P = 1 is spinning with 
an angular velocity uj = 20, embedded in an uniform background 
with p — P — 1. Again we used a glass like 3D particle distri- 
bution with 700 X 700 x 5 particles and periodic boundaries, and 
increased the amount of particles in the disk to achieve the desired 
density using the same mass per particle. For comparison, we ob- 
tained a reference simulation with ATHENA run at 400 x 400 cells. 
The results obtained with the Dedner scheme and with Athena 
are shown in Fig. (|3]). The shape, positions and amplitudes corre- 
spond very well, although the Gadget runs appears slightly more 
smoothed. Here, the cleaning scheme reduces the overall numerical 



3.2.3 Orszag-tang Vortex 

This planar test problem, introduced by lOrszag & TangI (Il979h , 
is well known for the interaction between several classes 
of shock waves (at different velocities) and the transition to 
MHD turbulence. Also, this test is commo nly used to vali- 
date MHD implementations (for example see iDai & WoodwardI 
19941: 'Picone & Dahlburd Il99ll: iLondrillo & Del Zaim^ l200d : 
Price & Monaghan 2005 : B0rve et albood) . It consists of ideal gas 
with 7 = 5/3 within a a box of (e.g. x = [0, l],^/ = [0, 1]) 
and periodic boundaries conditions. The velocity field is defined 
by Vx — —sin(2Tiy) and Vy — sin(27rx). The initial magnetic 
field is set to Bx = BqVx and By = Bosin{47vx). The initial den- 
sity is p = jP and the pressure is P = jBq. An usual time to 
evaluate the system t = 0.5. Fig. ^ shows the final result at that 
time for the magnetic pressure for the ATHENA run (left panel), 
and the Dedner scheme (right panel). The results are quite com- 
parable, however, the use of Spmhd method leaves its imprint in a 
slightly smoothed appearance in the Gadget results. This can also 
be seen in Fig. JS]), which shows a cut through the test for different 
implementations, comparable with other cuts done in the literature 
(B0rve et al. 2006). In general there is reasonable agreement, how- 
ever all the Spmhd results clearly show a smoothing of some fea- 
tures. However, the Dedner and standard implementations, tend to 
match better some regions that the dissipative schemes oversmooth 
(e.g. region near x 1 in Fig. [6]), even better that the Euler scheme. 
Note, that a exact comparison is difficult, mainly because this test 
includes the propagation of several types of magnetosonic waves, 
implying that if missing the correct velocity (i.e. by some dissi- 
pative effect) of a particular wave, the result will diverge between 
implementations. 

This periodic test in particular, is good for checking the 
Spmhd implementation based on the Euler potentials formal- 
ism, finding a very g ood agreement with other authors (i.e. 
iRosswog & Pricell2007h . As the Euler potentials are be V • B free 
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(a) Fast Rotor - Athena (b) Fast Rotor - Dedner 




Figure 4. Density distribution for the Strong Blast test at t = 0.02. The left panel shows the Athena results, while the right panel shows the Spmhd Dedner 
scheme. There is no visible difference between both plots, beside the small intrinsic smoothing of sharp features within the Spmhd solution. 



by construction, any numerical arising V B error can clearly traced 
back to the numerical inaccuracies in SPMHD formalism itself. Our 
major interest in this scheme is therefore the possibility to measure 
the errors that arise from the interpolation. 

In Fig. we show the calculated V • B errors as defined 
by Eq. It can be seen that the numerical errors in the standard 
Spmhd implementation are only slightly larger than the errors of 
the implementation based on Euler potential, however the spacial 
distribution varies. This is the numerical error limitation, which can 
be overcome by using higher resolution. Additionally, the numeri- 
cal V • B errors we see are caused by the magnetic field structures 
getting folded below the kernel scales. Then, the basic assumption 



on which Spmhd works, namely that the values of any quantity of 
interest are smooth below the kernel scales, starts to get violated. 
In Eulerian methods such structures are automatically mixed (e.g. 
dissipated) on the resolution scale. In Spmhd an extra scheme (i.e 
artificial dissipation) is needed to remove those small scale struc- 
tures, acting as regularization of the field below the kernel scale. 
The Dedner cleaning scheme acts in this way, however, as seen al- 
ready in previous test, it dissipates mainly the magnetic field struc- 
ture below the kernel scale and does not lead to strong smoothing of 
the field on scales larger than the kernel scale. However, is enough 
to significantly remove the numerical V • B errors, as shown in Fig. 
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Figure 6. A t = 0.5 cut of the pressure in the Orszag-Tang Vortex at 
y = 0.3125. The black line reflects the results obtained with Athena, 
while the colors show dif ferent MHD schem es. The cut can be compared 
with existing results, e.g. lB0rve et SI (l2006h . Note, that the variations in 
the solution between schemes mainly result from different dissipative char- 
acteristics. However, the overall solutions are in good agreement. 

More drastic approaches suffer from the same issue, as can be 
seen for the artificial dissipation case. 

3.2.4 Long time stability of the Orszag-tang Vortex 

The lo ng time evolution of the Or szag-Tang Vortex was first stud- 
ied by iPicone & Dahlburd (Il99lh . They focused on the possible 
stable solutions of the supersonic flow. Therefore, they studied the 
long time evolution varying initial Mach numbers. Interestingly, 
they always found a quasi stable configurations in the long time 
evolution, where they typically evolved the problem until t = 8. 
We run the Orszag-Tang vortex to large times to investigate the 
stability of the different implementations as well as the influence 
of the differently strong numerical dissipation. Additionally, we 
wanted to confirm th e limitations of the Euler potential formalism 
(Ib randenbur d |20 1 ol) . In Fig. [8] we show the long time evolution 
of the density distribution at various times. As expected, the im- 




Figure 7. V • B errors as defined by Eq. (|6) in the Orszag-Tang Vortex at 
t = 0.5. The upper row shows the standard scheme (left panel) and the 
Dedner scheme (upper right) panel. The lower row shows the calculated 
V • B errors from the Euler potential formulation (left panel) and a run 
with artificial dissipation (right panel). The Euler potential runs defines the 
real numerical limits of the simulation and confirms the need of a numerical 
cleaning scheme. 

plementation based on Euler potentials start to deviate from the 
expected solution quite early (t = 1) and even runs into some se- 
vere instability at larger times t > 4. The regularization scheme 
based on artificial dissipation, as well as periodically smoothing 
the magnetic field, show some significant effects of the underly- 
ing dissipation at times t > 4. These deviations even develop an 
instability within the scheme based on artificial dissipation, and 
even earlier times t < 4 for higher dissipation constants. How- 
ever, both the standard MHD implementation as well as the one 
based on the Dedner cleaning scheme show excellent performance 
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in the long term e volution and stability, well comparing to the re- 
sults presented in IPicone & Dahlburgj (Il99lh . It is worth to note 
that both schemes appear to have less numerical dissipation than 
Athena (See http://www.astro.princeton.edu/'-^jstone/athena.html 
for a comparison run until T = 5). In our ATHENA run, the two 
central density peaks start to approach each other, while they are 
still stable in both, the standard as well as the Dedner Spmhd im- 
plementation. 

3.3 Discussion 

In Fig. ([9]) we show the comparison of the performance of the dif- 
ferent Spmhd implementations among all the tests discussed in 
the previous sections. We compare the numerical errors obtained 
by the different implementations against the standard implementa- 
tion. Additionally, we re- size the points, follow ing the quality cri- 
terion as defined in IPolag & Stasvszvnl (l2009h . Therefore we are 
able to compare the performance between all against the standard 
runs. For this, we define 3 different sizes when comparing to the 
standard quality value, corresponding to different grades in the im- 
provement. 

The Dedner cleaning schemes clearly reduced the V • B er- 
ror in all tests presented. Although the dissipative schemes have 
even lower V • B errors, the cleaning scheme out stands, as it low- 
ers the V • B and also is least dissipative. In contrast, there some 
tests showing lower V • B errors, but the regularization schemes 
over smooth some features, therefore enhancing differences with 
the correct solution. 

Especially regularization schemes like the one based on arti- 
ficial dissipation can be seen to mimic the Ohmic dissipation and 
therefore leading away from ideal MHD. In galaxy clusters, such 
effects might have to be taken into account to obtain results which 
better agree with the observed profile of the resulting magnetic 
fields teonafedee t al. 2011), and therefore it is quite important to 
have an underlying MHD implementation which does not suffer 
from any form of artificial dissipation. To study the structure of 
the magnetic field as imprinted by the complex, hydro-dynamical 
flows as imprinted by structure formation, it is quite important to 
have a scheme which on one hand regularizes the magnetic field be- 
low the kernel scale to avoid unwanted numerical artifacts as well 
as not influencing the magnetic field structure at scales larger than 
the kernel scale. As shown by the tests performed in the previous 
section, the Spmhd implementation based on the Dedner V • B 
cleaning scheme seems to fulfill these requirements. 



4 GALAXY CLUSTER AND MAGNETIC FIELDS 

In the hierarchical picture of structure formation, small objects 
collapse first and then merge to subsequently form larger struc- 
tures. This formation process is reflected in the intricate structure 
of galaxy clusters, which properties depend on how the structure 
of the smaller objects was accreted by the cluster and its gravita- 
tional potential. These accretion or merging events cause shocks 
and turbulence in the ICM, leading to a redistribution or amplifi- 
cation of magnetic fields. To describe this process within numer- 
ical simulations is extremely challenging, as the structures in and 
around galaxy clusters are intricate and range over many dynam- 
ical scales and orders of magnitude. Here, we used the so called 
Zoomed Initial Con ditions (ZIC) technique (iKatz & Whit3ll993l : 
Tormen et al ][i993) to follow the evolution of a single, relatively 
low mass galaxy cluster within a large cosmological box with so far 



unreached precision to study the evolution and the final structure of 
the magnetic field and compare it with the observed magnetic field 
structure within a similar, low mass galaxy cluster. 

To perform this study, we re-simulated a Lagrangian region se- 
lected from a lower resolution dark matter only cosmological box. 
This parent simulation has a box-size of 684 Mpc, and assumed a 
flat ACDM cosmology with Vtrn = 0.3 (matter density), h = 0.7 
(Hubble constant), fbar = 0.13 (baryon fraction) and as = 0.9 
(normalization of the power spectrum). The selected cluster with 
a final mass of 1.5 x 10^^ Mq was re-simulated using 4 different 
particle masses for the high resolution region. To optimize the setup 
of the initial conditions, the high resolution region was sampled 
with a grid of 16^ cells, where only sub-cells were re-sampled at 
high resolution to allow for quasi arbitrary shapes. The exact shape 
of each high-resolution region was iterated by repeatedly running 
dark-matter only simulations, until the targeted objects were clean 
of any lower-resolution boundary particle out to 3 5 virial radii. 
The initial particle distributions, before adding any Zeldovi ch dis- 
place ment, were taken from a relaxed glass configuration (IWhitd 
Il996h . The highest-resolution (e.g. 130x) thereby corresponds to 
a mass of the dark matter and gas particles of 1.2 x lO^M© and 
2.3 X 10^ M© respectively. The according gravitational softening 
is 1.6 kpc. Table m summarizes the settings for the different reso- 
lutions. For simplicity we assumed an initially homogeneous mag- 
netic field seed of 10~^^ G. Additionally, we ran several of the low- 
resolution versions of the simulations with our different SPMHD 
schemes to check the influence of the different techniques on the 
results. To focus on the shaping of the magnetic field structure due 
to the structure formation process, we ran the simulation without 
including any additional physics like cooling or star-formation. 

4.1 Galaxy duster slices and profiles 

To show the complexity of the cluster atmosphere we cut thin slices 
through the center of the highest resolution version of the simula- 
tion (130x), centered on the potential minimum within the galaxy 
cluster at 2; = in Fig. ( [Tol ). Clearly, the accretion shock at the 
outskirts (r > 2 Mpc) is visible, getting penetrated by cold, fila- 
mentary structures plunging into the hot, cluster atmosphere. Also 
several internal shocks (from previous merger events) are still vis- 
ible within the outer parts of the atmosphere, just approaching the 
virial radius. In contrary, the central part of the cluster looks rela- 
tively relaxed. In a high /3 plasma as the simulated here, the mag- 
netic field lines are frozen into the velocity field, while being com- 
pressed and advected together. Therefore, especially beyond the 
virial radius, the structures seen in the magnetic field are well cor- 
related to the slightly denser gas, left over material of several thou- 
sands of resolved substructures, where the gas was stripped while 
the sub- structures got dissolved within the denser environment of 
the cluster during its assembly. Within the virial radius, structures 
are more virialized and mixed. Additionally, filamentary structures 
within the magnetic field can be noticed due to the presence of tur- 
bulence. Also, additional amplification at the shock fronts, with a 
mild decline afterwards as expected, are clearly visible. The distri- 
bution of the numerical V • B errors shows that they are mostly 
located in the very low density region, where the resolution of the 
underlying particle distribution is quite low. 

In Fig. (fni we show the calculated radial profiles of the mag- 
netic field. With enough resolution (to resolve the turbulence prop- 
erly), the magnetic field saturates at roughly equipartition of an 
fraction of the turbulent energy fraction (a fraction of 10 percent 
would lead to values of order of lO/xG). Naturally from thermody- 
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Figure 8. Density maps of the Orszag-Tang vortex at various times T = 0.5, 1,2, 4, 6 using several MHD schemes. There are deviations from the expected 
solution (Athena), even appearing at early times T = 1.0 for Euler, and later in the more dissipative schemes. There are also problems long-time simulations. 
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Table 1. Overview of the different initial conditions. The last column shows the final total number of particles particles inside the virial radius. 



namics, the energy content of the magnetic field should be com- 
parable to the turbulent energy, which is also observed in galax- 
ies. This is the result of a turbulent dynamo operating, exten- 
sively demonstrated by numerical simulations of colliding galaxies 
jKotarba e t al. 2011; Geng et al. 201 2) and simulation s of galactic 
halo formation (iBeck et aDl2012h . In lBeck et aD (1201 2h it was even 
shown, that the evolution of the magnetic field observed in simula- 
tions very well can be described by a simple, turbulent dynamo 
model. However, in galaxy clusters the observed magnetic field 
is found only to be order of /xG, so significantly below equipar- 
tition values from the level of turbulence present in galaxy clusters. 



Consistent with previous studies jBonafede et al.ll201ll) we found 
that without explicitly modeling physical dissipative processes, the 
magnetic field saturates at larger values than indicated by the ob- 
servations. It is important to stress here, that the numerical V • B 
errors are decreasing with increasing resolution (Fig. dTT}), while 
the magnetic field profile stays converged. Also, with the exception 
of the lowest resolution run (Ix), the results of the standard Spmhd 
simulations compare well with the Dedner implementation, regard- 
less of the fact that the Dedner shows less numerical V • B error. 
In general, within the 130x simulation, the numerical V • B error is 
even more than one order of magnitude lower than for the 6x run. 
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Figure 9. Summary plot showing the results of all ID and 2D tests for the different Spmhd implementations, comparing the average | V • B| errors from the 
different schemes (symbols) and the Dedner V • B cleaning against the st andard Spmhd impl ementation. In black are the shock tube tests and in colors the 
advanced problems. The points sizes follows a quality criterion defined in'Dolag & Stasvszvnl (2009), comparing the Spmhd solutions with the Euler ones. 
We define 3 different sizes when compared to the standard quality value, corresponding to worst (small), major (big) or minor (normal) improvement. 



while leaving the radial magnetic field profile unchanged. This is 
another strong prove that the results are not driven by numerical 
artifacts. 



4.2 Synthetic RM Measurements 

To compare with observations and especially to study the predicted 
magnetic field structures with in the simulations , we produced syn- 
thetic RM maps, using Smac jPolag et al.l2005h , which performs a 
line of sight integration and projection of the simulation SPH data 
onto a grid. The simulated cluster was not constrained to repro- 
duce a certain, observed cluster in particular and we in any case 
neglect the effect of a physically motivated magnetic dissipation, 
which would be needed to match the exact radial magnetic profile 
of galaxy clusters. The simulations also neglect the effects of cool- 
ing, star-formation and stellar, as well as AGN feedback, which 
would be needed to reproduce certain aspects of the thermal struc- 
ture within galaxy clusters. Therefore, we do not expect to match 
the exact amplitude of the RM signal of the cluster and therefore 
investigate only the structural properties of the obtained RM maps. 

Fig. iT2\ gives an impression of the RM maps obtained from 
the Ix resolution simulations. The structural properties of the RM 
patterns agree very well between the standard implementation and 
the implementation based on the Dedner scheme. We here show 
also the results obtained from the Euler Potential implementation, 
which passively (e.g. neglecting the Lorenz force) evolves the mag- 
netic field in the simulation. As this implementation just reflects the 
integrated winding of the magnetic field by the hydrodynamic pat- 
terns within the ICM, it demonstrates that the magnetic field pat- 
terns observed in galaxy clusters is strongly related to the underly- 
ing turbulence within the ICM. We also show the effect of numer- 
ical dissipation and magnetic regularization on the predicted mag- 
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Figure 13. Synthetic RM maps for different resolutions with the Dedner 
scheme. The maps show the central region (174 x 174 kpc) of the simulated 
cluster. From upper left to lower right: Ix, 6x, lOx and 130x resolution. 

netic field structures, it can be seen that numerical magnetic dis- 
sipation can lead to severe smoothing of the magnetic field within 
such simulations. 

In Fig. ( [T3] ) we show synthetic maps obtained with the Dedner 
implementation at the different resolutions. With increasing reso- 
lution, more magnetic field reversals are resolved. This is expected 
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Figure 10. Slices through the center of the galaxy cluster at redshift 2; = for 130x resolution, showing the large scale cluster environment. Shown are the 
density, temperature, magnetic energy and the V • B errors. The dashed circle indicates the virial radius of the cluster. 
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Figure 11. Radially averaged, magnetic field profiles at z = for different resolutions (left panel) obtained with the Dedner implementation. The colored 
lines correspond to different resolutions. Numerical profiles of the V • B errors (right panel). Additionally we show the 130x run without cleaning scheme in 
dotted lines. Increasing resolution yields lower errors and a better resolved magnetic field. 
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Figure 12. Synthetic RM maps of the central region of galaxy clusters from the Ix run, using various MHD implementations. From left to right: standard 
Spmhd, Dedner cleaning, artificial dissipation, B smoothing and Euler potentials without Lorentz force. Shown is a projected region of 174 x 174kpc centered 
onto the cluster. 



as with higher resolution we are able to resolve smaller structures 
within the turbulent velocity field, which twist and bend the mag- 
netic field. Note, that with the highest resolution (e.g. the 130x) we 
reach a resolution, which comes close to the coherence lengths of 
the observed magnetic field on clusters of a few kpc. For exam- 
ple, the radio lobes of 3C449 (iFeretti et al.lll995h can be used to 
probe the magnetic field of the cluster atmosphere by its RM im- 
print. The radio source 3C449 is located in the center of a relatively 
low mass galaxy cluster with a observed temperature of the ICM of 
Ikev (Hardcastle et al. 1998|) at z = 0.1711. Fig. compares 
the observed RM patterns with the RM patterns obtained from the 
simulations, cutting the synthetic maps to the shape of the observed 
radio lobe, mimicking the observational window. Both, simulations 
and observations not just show random patterns but also some more 
filamentary structures as expected from the presence of MHD tur- 
bulence. 

4.3 Structure Functions 

To take the observational limitations into account when compar- 
ing the results of our simulations with observations, we used syn- 
thetic observations and masked them with the shape of the Faraday 
screen, as done in Fig. ([T4l ). This mimics any kind of bias intro- 
duced by the observational window - defined by the size and shape 
of the observed radio lopes - for our synthetic calculations. The ob- 
served RM maps also sho\ y a significant cons tant offset caused by 
the galactic magnetic field iFeretti et al.l (Il995h , so we subtract this 
galactic contribution from the observations when comparing with 
theoretical predictions. 

To obtain a quantitative comparison we calculated the pro- 
jected second order structure function 

S(t) = (^(RM(rO -RM(/ + r))^^, (7) 

of the observed and synthetic RM maps, with r = yAx^+A^ 
being the distance from a pixel at position r' — (x,y). The re- 
sulting matrix is then averaged in radial bins to obtain the structure 
func tions. 

lKraichnanl ( ll965h formulated one of the first phenomenologi- 
cal theory of MHD turbulence. He found that in the presence of a 
strong mean magnetic field, will result in an magnetic energy spec- 
trum of ^ k~^^^. However, if the magnetic field is weak, as used 
in most of the mean field theories, the expected magnetic energy 
spectrum should be Kolmogorov like and ^ k~^^^. Both theo- 
ries assume isotropic turbulence, which might not the case in as- 
trophysical systems, leading to a suppression of the energy cascade 



Obs 




Figure 14. Comparison of the predicted RM map from simulations at 
the highest resolu tion (e.g. 130x) with the observed map of 3C449 in 
IFeretti et al] ( Il995h . The synthetic map is cut to the same shape as the ob- 
servational map to allow a better comparison. Both maps have a resolution 
of approximately 0.17kpc per pixel. 

along the direction of the mean magnetic field. However, observa- 
tions on galactic sca les jHan et al. 2004) and numerical simulations 
jMason et al.ll2008h . favor the Kolmogorov spectrum and energy 

transfer from l arger to small scales. 

Following lRuzmaikin. A. A. and Sokolov, D. D. and Shukurov A.M.I 

('1988!), rough estimates for the structure functions for the RMs can 
be derived. There is a coherence length /o, determining when the 
RM are not longer specially correlated. Thus, for scales d ^ lo 

S{d)c^RMl(^^y (8) 

with values of 7 = 2/3 for the Kolmogorov spectrum and 7 = 1/3 
for the Kraichnan spectrum. At scales larger than Iq the structure 
functions should have a constant value RMq and the RMs should 
be uncorrected. 
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Figure 15. Structure function derived from the observed rotation mea- 
surement (black line) of 3C449 in Feretti et al. (1995) compared to the 
expectation of simple, theoretical models. In red the expected shape for 
the Kolmogorov spectrum and in blue the one expected from a Kraichnan 
spectrum. The dotted lines correspond the parameter values of RMq = 
35 rad and lo = 20 kpc, whereas the dashed lines correspond to 
RMo = 38 rad and Iq = 45kpc. Note, that different values for the 
parameters allow for different fits of the slopes. 



Figure 16. Structure functions calculated from the observed RM map (black 
line) of 3C449 in Feretti et al. ( 1995) compared to 130x simulations using 
the Dedner implementation (red lines). The dashed line shows the structure 
function calculated from the full (174 x 174 kpc) maps (see Fig. (fT3l), 
whereas the solid lines show maps cut to the shape of the observed RM 
map (as shown in Fig. (flU). The dotted Hnes indicate changes in the slope 
at around a scale of 10 ~ 20kpc, visible in the structure function obtained 
from the simulations. 



In Fig. ( [T5] ) we compare theoretical predictions for the struc- 
ture functions of those simplified models with the structure functio n 
computed from the observed RM of 3C449 iFeretti et al] (Il995h . 
In red is shown the expected shape for a Kolmogorov like spec- 
trum and in blue the expectation for a Kraichnan spectrum. Un- 
fortunately, the current observations do not allow to distinguish 
these two models. The Kolmogorov like spectrum seems to fit 
slightly better, in line with the findings for the ma gnetic field power 
spect rum reconstructed from the RMs in Hydra (iKuchar & EnBlinI 
l201lh . But at larger distances the RM should be uncorrelated and 
therefore have a constant structure function. As the spatial range 
over which the magnetic field can be measured is still limited to 
the core region of the cluster, density and magnetic field are not 
significantly declining over the region covered by the radio lobes 
and therefore can not explain the decrease of the structure func- 
tion at larger distances. However, the synthetic RM maps and the 
windowing introduced by the shape of the RM region where data 
is available, can responsible for the radial decline of the structure 
function at these distances (see Fig. [16]). 

A quantitative comparison between the structural properties 
of the observed RM and structures obtained from the simulations 
is shown in Fig. ([T6l) . As mentioned and discussed earlier, we were 
interested in the shape of the structure function and therefore re- 
normalize the amplitude of synthetic RM to the same, global RM 
value of the observations to make them comparable. The shape of 
the structure function at the highest resolution is relatively flat at 
scales above 10 ^ 20kpc, and shows a steepening at scales smaller 
than that. The extreme resolution of our simulation is reflected in 
the fact, that the structure function does not drop to zero even at the 
smallest scales showing still significant structures. Especially when 
the calculation is restricted to the spatial observed mask in 3C449 
it matches the observations very well. How the results change for 
different diffusive numerical schemes or at different resolution is 
shown in appendix [A] As soon as simulations are reaching reso- 
lutions similar to the observed scales, the magnetic field structure 
obtained by the simulated galaxy clusters, explains well the mag- 
netic field structures observed through RMs. This means that struc- 
ture formation, by its hierarchical nature and the related turbulence 



within galaxy clusters is already shaping the magnetic field struc- 
ture within galaxy clusters to observed structures. 



5 CONCLUSIONS 

We continue the development started bv lPolag & Stasvszvnl (l2009h 
of an Spmhd implementation of MHD in the cosmological simula- 
tion code Gadget. We performed various standard test problems 
and discussed instability corrections, regularization schemes and 
V • B cleaning techniques. The code was applied to cosmologi- 
cal astrophysical simulations, studying the role of resolution and 
regularization schemes in simulations of galaxy clusters. The main 
focus was set on the role of the V • B = constrain, also compar- 
ing with Eulerian solutions. Our main findings can be summarized 
as follows: 

• Correcting the tensile instability by explicitly subtracting the 
contribution of a numerical non-zero divergence of the magnetic 
field to the Lorenz f orc e from the Maxwell tenso r as suggested 
by B 0rve et al.l (1200 ih or lPolag & Stasvszvnl (l2009h is performing 
well. To avoid spurious effects due to sampling problems, particu- 
larly in the front shocks, a threshold in the correction improves the 
performance and quality. This feature turns out to be fundamental 
for simulations with extremely high spatial resolution, where pre- 
viously the numerical instabilities dominated and the simulations 
could not have been performed. 

• We successfully implemented a multidimensional divergence 
cleaning method in Spmhd. The GADGET Spmhd implementa- 
tion continues to perform very well on our multidimensional shock 
tube tests as well as on commonly used planar test problems. We 
showed that the Dedner cleaning scheme does not affect the shape 
of shocks, but reduces the V • B errors. This is important in as- 
trophysical situations. Also, the evolution of the cleaning is com- 
pletely local and only affected by the Spmhd interpolants. 

• Testing our schemes in long time runs of the Orszag-Tang Vor- 
tex, we found that regularization schemes which depend on dissipa- 
tion as well as implementations based on Euler potentials can lead 
to instabilities, whereas the standard implementation as well as the 
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Dedner cleaning scheme are found to be robust and stable. This is 
a warning on the limits of possible MHD implementations within 
astrophysical environments. 

We obtained our most important findings by comparing with 
an Euler (divergence-free by construction) scheme. This shows, 
that we have already reached the V • B error levels inherent from the 
Spmhd scheme itself (see Fig. [7]), and any further cleaning or reg- 
ularization will require to numerically dissipate the magnetic field. 
However, this does not yield an error-free implementation, thus still 
some V • B related instabilities can arise. In the case of the Euler 
potentials, we stress the fact that this representation lacks important 
features in the evolution of the magnetic field (Brandenburg 2010). 
However, outstanding progress has been made on providing stable 
MHD and Spmhd implementations, removing unwanted instabili- 
ties. 

Furthermore, our code has different features with each reg- 
ularization or cleaning scheme used. The success of the Dedner 
cleaning brings the simulations closer to the "ideal" MHD state, by 
locally calculating and subtracting the error term in the induction 
equation. 

In Fig. ^ we show a comparison between all the implemen- 
tations in different tests. The dissipative schemes have lower V • B 
errors, however the Dedner cleaning schemes stands out, lowering 
the V • B error, and can in most cases improve the results. In con- 
trast, there are some tests showing that to lowering V • B errors can 
over-smooth important features in contrast to a correct solution. 

We ran cosmological simulations of galaxy clusters with dif- 
ferent underlying SPMHD schemes for the magnetic field evolu- 
tion. The magnetic field profiles of the different simulation schemes 
show similar properties, even when comparing with Eulerian codes. 

The comparison between the different schemes leads us to the 
conclusion that in general, the non- vanishing numerical V • B terms 
are not important for the proper evolution or growth of the magnetic 
field. 

It is known, that in galaxy clusters, the growth of the 
magnetic field is not only driven by the adiabatic col- 
lapse. The magnetic fields lines are strongly tangled on 
all scales, leading to a complex, turbulent morpholog y 
(iRuzmaikin. A. A. and Sokolov. D. D. and Shukurov A.M.! 1 19881) . 
An intrinsic feature of Spmhd is its adaptivity, particularly in 
spatial resolution. Therefore, the turbulent cascade is well de- 
scribed also within the framework of the collapse of cosmological 
structures, leading to a physical growth of the magnetic field by a 
turbulent dynamo. Note, that previous work clearly demonstrated 
that equipartition between the magnetic energy and the turbulent 
energy is reached, independent of the magnetic seed field strength. 
Therefore, the predicted magnetic field structure within galaxy 
clusters can be studied independently of the real origin of the 
magnetic seed fields. By analyzing synthetic RM maps, we 
evaluated that previously used dissipative regularization schemes 
are producing numerically driven, large reversal scales within 
the magnetic field, which are not comparable to what is seen in 
observations. However, our standard Spmhd implementation as 
well as our Dedner implementation allow to follow the magnetic 
field structures without loosing resolution by dissipative effects. 
In the extremely high resolution simulations we performed, 
magnetic field structures down to kpc scales are for the first time 
resolved within cosmological simulations of structure formation. 
By studying structure functions of synthetic RM maps, we found 
that these simulations reproduce very well the structures within 
observed RM maps. This demonstrates that the observed structure 



of the magnetic field within galaxy clusters can be shaped already 
by the cosmological structure formation process and the thereby 
driven, turbulent gas motions. 

Finally, the use of advanced Spmhd schemes allows to simu- 
late the properties of the ICM at very high spatial resolutions, con- 
tributing to the global picture of the properties of the hot gas within 
galaxy clusters. 
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APPENDIX A: DEPENDENCE OF THE STRUCTURE 
FUNCTION ON THE NUMERICAL SETUP 

Comparing the structure functions obtained from simulations using 
different numerical implementations and underlying resolutions, 
demonstrates how effective these simulations can describe the de- 
tailed magnetic field structure observed in the RM maps of galaxy 
clusters. 

In the case of comparing different implementations, we can 
study the effective resolution reached due to the numerical dissipa- 
tion within the different implementations. 

For each simulated cluster we calculated the synthetic RM 
maps with a physical size of 174 x 174 kpc (see Fig. [12]) and 
calculate their structure function shown in Fig. dAlb , ignoring 
the windowing which would be introduced by the visibility of the 
real observations. To account for the varying efficiency to resolve 
the underlaying dynamo mechanism among the different numerical 
implementations and to make them comparable with the observa- 
tions, we rescaled the synthetic maps using the as a reference the 
amplitude at large scales in observations, assuming that at those 
scales the RMs are uncorrected. For the schemes, where the reg- 
ularization is based on dissipation, the effective resolution is dras- 
tically reduced by large factors compared to the SPH smoothing 
scale. Therefore the resolved coherence length at the Ix resolution 
is reduced to scales comparable to the size of the maps used and 
manifests itself in a structure function which still keeps growing. 

In Fig. (IA2I) a resolution study using the Dedner implemen- 
tation is shown by the structure functions calculated from the 
174 X 174 kpc RM maps. As before, we scaled the synthetic maps 
using the amplitude of the observations. Increasing the resolution 
we are able to resolve smaller scales in the underlying turbulence, 
reflected in being able to reach smaller and smaller scales within the 
RM maps. Also, the coherence length of the turbulent field shifts 
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whereas the small scales are not affected. Here, the simulations us- 
ing the Dedner cleaning, at the highest resolution, show a remark- 
able match to the observations. However, the simulation without 
the Dedner cleaning scheme, shows a clear excess of small scale 
structures. The variability at small scales, indicate the need to in- 
clude additional physics (i.e. ohmic dissipation, turbulent diffusion 
or viscosity) when further increasing the resolution. This will be 
needed to limit the turbulent cascade to not extend below the scales 
where the fluctuations are present in the observations. 



Figure A2. Shown are structure functions S calculated from the full (174 x 
174 kpc) maps (see Fig. [13) for the simulations with increasing resolution 
using the Dedner implementation and in dashed the highest resolution for 
the standard case (i.e. without cleaning scheme). As a reference we show 
the result of the from the observed rotation measurement (doted black line). 
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Figure A3. Shown are structure functions S cal culated from the ob served 
rotation measurement (black line) of 3C449 in Feretti et aP (Il995h com- 
pared with the simulations at different resolutions with the cleaning scheme. 
The structure functions are calculated from the map which has been cut to 
the shape of the region of the observed rotation measurements (as shown in 
Fig. [14). The dashed blue line shows the 130x resolution run, without any 
cleaning scheme. 



to smaller and smaller scales, as visible by the beginning of the 
flat part moving towards smaller scales when increasing resolution. 
There are no considerable differences between the 6x and lOx runs, 
as the increase in resolution is quite minor in that case. This demon- 
strates the need to vary the resolution by a large amount when in- 
vestigating resolution effects, as done here by increasing the reso- 
lution overall by a factor of 130. In the 130x case, we also show 
the case without Dedner cleaning, where the power at small scales 
is even larger and the coherence scale is moved very close to the 
resolution scale. Here it can not be excluded, that these fluctuations 
within the magnetic field could be driven by the residual numerical 
V • B errors within the simulations, which do not make use of the 
Dedner cleaning scheme. 

Finally, Fig. (|A3j, shows the result where we calculated the 
structure function for the maps, cropped to the shape of the ob- 
served rotation measurements (as shown in Fig. [14). The win- 
dowing due to the shape of the actual observations leads to a de- 
cline of the correlation function at large scales, bringing the simu- 
lated structure functions in good agreement with the observed one. 
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